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Abstract 

Gravitational lensing of a quasar by a spiral galaxy should often be accompanied by 

Z^ ' damped Lya absorption and dust extinction due to the intervening gaseous disk. In nearly 

edge-on configurations, the surface mass density of the gas and stars in the disk could by 

^\1 ' itself split the quasar image and contribute significantly to the overall lensing cross section. 

We calculate the lensing probability of a disk+halo mass model for spiral galaxies, including 

C*^ . cosmic evolution of the lens parameters. A considerable fraction of the lens systems contains 

0^ ' two images with sub-arcsecond separation, straddling a nearly edge-on disk. Because of 

(-H ^ that, extinction by dust together with observational selection effects (involving a minimum 

O , . separation and a maximum fiux ratio for the lensed images), suppress the detection efficiency 

of spiral lenses in optical wavebands by at least an order of magnitude. The missing lenses 

could be recovered in radio surveys. In modifying the statistics of damped Lya absorbers, 

the effect of extinction dominates over the magnification bias due to lensing. 



1 Introduction 



Gravitational lensing by a spiral galaxy occurs when the line-of-sight to a background quasar 
cH ' passes within a few kpc from the center of the galactic disk. Since galactic disks are rich in neutral 

hydrogen (HI), the quasar spectrum is likely to show a damped Lya absorption trough at the lens 
redshift. Therefore, the efficiency of blind searches for gravitational lensing with sub-arcsecond 
splitting can be enhanced by more than an order of magnitude, by selecting a subset of all bright 
quasars which show a low-redshift {z '^ 1) damped Lya absorption with a high HI column density, 
> 10^^ cm^^ (Bartelmann & Loeb 1996). Moreover, multiply imaged quasars could be identified 
spectroscopically through their multiple-trough absorption spectrum. The composite spectrum of 
a lensed quasar is a superposition of the spectra received from the different images which intersect 
the absorbing disk at different locations, probe different HI column densities, and hence acquire 
different widths in their damped Lya troughs (Loeb 1997). 

The magnification bias due to lensing changes the statistics of damped Lya absorbers (DLAs) 
in quasar spectra by bringing into view quasars that would otherwise fall below the detection 
threshold (Bartelmann & Loeb 1996; Smette, Claeskens, & Surdej 1997). For optical observations, 
this effect is counteracted by dust extinction in the lensing galaxy (Malhotra, Rhoads, & Turner 
1997; Perna, Loeb, & Bartelmann 1997). The combination of lensing and dust extinction results 
in a net distortion of the HI column density distribution of damped Lya absorbers. 

Since galactic disks are thin, their own surface mass density could exceed the critical value 
necessary for image splitting when they happen to be projected nearly edge-on. In such orienta- 
tions, a pair of quasar images will straddle the symmetry axis of the disk rather than its center. 
When averaged over all possible disk orientations, the total lensing cross section of the disk-|-halo 



mass distribution should still be roughly the same as that of a halo with the same spherically- 
averaged mass profile (Wang & Turner 1997; Keeton & Kochanek 1997). However, because of the 
new predominant image configurations around edge-on disks, the effects of dust extinction and HI 
absorption are substantially different when Icnsing by the disk is included in the calculation. 

Previous theoretical studies have either ignored the possibility of lensing by the disk itself 
(Bartelmann & Loeb 1996; Smette et al. 1997), or else ignored the effects of dust extinction or HI 
absorption on the lensing statistics (Mailer, Flores, & Primack 1997; Wang & Turner 1997; Keeton 
& Kochanek 1997). In addition, these discussions did not examine the implications of evolution in 
the spiral lens population on its overall lensing properties and statistics. In this paper, we include 
all of the above ingredients. Our model evolves the number density of galaxies based on the Press- 
Schechter (1974) theory for dark matter halos, combined with simple evolutionary scaling laws for 
the disk properties inside these halos (Mo, Mao, & White 1997). 

In §2, we describe our model for the mass and HI distributions in galactic disks, as well as our 
model for their dust content. We also summarize the adopted scaling laws for the evolution of the 
lens population with redshift. In §3, we compute the lensing cross sections, imaging probabilities, 
and lensing statistics of our model spiral galaxies. The impact of dust on the lensing statistics, 
and the effect of lensing on the statistics of damped Lya absorbers, are examined in §4. Finally, 
§5 summarizes the main conclusions from this work. 

2 A Lens Model for Spiral Galaxies 

2.1 Mass Model and Lensing Properties 

We use a mass model for spiral lenses which was recently suggested by Keeton & Kochanek 
(1997). The model consists of oblate, ellipsoidal, isothermal building blocks for individual mass 
components. They have the axisymmetric density distribution 

^^^' '^ ^ 4^Gg3 sin-i e i?2 + ,2 + ^2/^2 > (1) 

where R is the distance from the symmetry axis of the ellipsoid, 2:3 is the distance from its mid- 
plane, (73 is the axis ratio, s is the core radius which softens the mass distribution, Vc is the 
asymptotic circular velocity, and e is the eccentricity of the mass distribution, 

e - (1 - cilf'^ ■ (2) 

Hence, each mass component is described by one parameter quantifying its dynamical properties, 
namely the circular velocity Uc, and two parameters characterizing its shape, namely the core 
radius s and the oblateness (73. 

The projection of the three-dimensional density p results in the two-dimensional surface mass 
density 

nx)^^h,[q\s' + xl)+xl]-"\ (3) 

Here, Scr is the critical surface mass density for lensing, 

c^ D 

(4) 
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where D\, Ds, and D\s are the angular-diameter distances between the observer and the lens, the 
observer and the source, and the lens and the source, respectively, x — (a;i,a;2) is the position 
vector on the sky, and 

&p = 27r( — 1 — r— 1- =^ . -1 , (5) 

\ c / Us sm e sm e 

where h is the Einstein radius of a singular isothermal sphere with a circular velocity Wc. Finally, 

q is the projected axis ratio 

q={qlcos^i + sin^if'^, (6) 



where i is the incHnation angle of the disk relative to the line-of-sight, with i = for an edge-on 
orientation. 

As shown by Keeton & Kochanek (1997), the lensing potential of the surface- mass density (H) 
is 

„2ll/2 



■0(5, 93) = xiai + X2a2 - bpS In [{g + s)^ + (1 - q^)xl] 



(7) 
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r 2/ 2 I 2\ I 2 

= 1.9 {xi +s )+ .T2 



21 1/2 



a2 



(1-92)1/2 

(l-g2)l/2 



tan 



tanh 



(1 



g2)i/2^i 



g + s 

(l-g2)l/2 



X2 



Q + q^s 



(8) 



We now combine three of the isothermal oblate ellipsoids to a mass model consisting of a 
maximal disk and a surrounding halo. The total lensing potential then reads 



f/' = V'(sd, gsd) - '0(''d, gsd) + V'(sh, 1) 



(9) 



where Sd,h are the core radii of disk and halo, respectively, q^ is the disk axis ratio, and rd is the 
disk truncation radius, or disk radius for simplicity. The first term is the potential of a disk with 
an asymptotically flat rotation curve, axis ratio gad, and core radius s^. The second term truncates 
that disk at radius r^. The third term adds the surrounding spherical halo necessary to maintain a 
flat rotation curve beyond the disk truncation radius. This is the maximal truncated Mestel disU^ 
model introduced by Keeton & Kochanek (1997). All lensing properties of the combined model 
can now be calculated in terms of the potential ip. 




Figure 1 : Rotation curve of a model consisting of a truncated Mestel disk embedded in an isother- 
mal halo with Vc = 220 km s"'^. The curve steeply rises to v^ and then levels off, being flat to 
better than 5% at large radii. 



To completely define the model, we need to specify five parameters, namely the circular velocity 
Vc, the two core radii Sd,h, the disk radius r^, and the disk oblateness q^d- For reference, we define 
Vc* — 220kms~^, r^* — 8/i~^kpc, and q^* = 0.03. Requiring a flat rotation curve, we must then 
use Sh* ~ 0.72 rd* (Keeton & Kochanek 1997). Finally, we choose Sd* — 0.2 h~^ kpc. The rotation 
curve of this model is plotted in Figure |l|. As the figure shows, the rotation curve is flat to better 
than 5% beyond a radius of ~ 2 h~^ kpc. 

To illustrate the lensing properties of the combined disk-|-halo model described by equation (|9|) , 
we show in Figure g the caustics and critical curve configuration for a disk inclination angle of 
i = 10°. The lens has two critical curves and caustics. The inner caustic, which is the image of 
the outer critical curve, has the familiar astroid shape with four cusps connected by folds. 

iPor S(i — > 0, the disk becomes a Mestel (1963) disk. 




Figure 2: Critical curves (dashed lines) and caustics (solid lines) for a disk+halo lens. Since near 
edge-on configurations are preferred both geometrically and by the magnification bias, we show 
results for a disk inclination angle of i = 10°. The caustic corresponding to the outer critical curve 
shows the familiar astroid shape with four cusps connected by folds. The substantial contribution 
of the disk to the surface mass density stretches the caustic along the symmetry axis of the 
projected disk. 

2.2 Neutral Hydrogen Disk 

Next, we address the effects of the gaseous component of the disk. To evaluate the level of dust 
extinction and HI absorption, we assume that the HI density, nn, follows a double-exponential 
profile, 

\^3\\ ( R 



nii{R,X3) = nH,o exp 



(10) 



The three parameters that define this profile are the central neutral hydrogen density hh.Oi ttic disk 
scale height H, and the disk scale length Rq. In order to keep the number of free parameters small 
in our model, we assume H = q^r^ and Rq = rd, so that the HI scale height is given by the short 
axis of the oblate ellipsoid used to describe the lensing disk. We also identify the HI scale radius 
with the truncation radius for the lensing disk. The radial shape of the resulting radial HI profile 
is similar to that observed in local disk galaxies (Broeils & van Woerden 1994). We choose nH,o 
so that the face-on HI column density, A^o — "^Hnufi, corresponds to the characteristic observed 
value (Broeils & van Woerden 1994), Nq « 11.25 x 10^°cm~^. Hereafter, we use the notation 
N = 10^° A^20 cm~^. Note that with these parameter choices, the solar neighborhood values for 
the HI column density (A'^hi = 7.4 x 10^° cm"^; Kulkarni & Heiles 1987) and HI scale height (a 
few hundred pc; Knapp 1987; Kuijken & Gilmore 1989) of the Milky Way disk are reproduced 
reasonably well. 



2.3 Extinction by Dust 

We assume that the distribution of dust follows equation (noh with a constant ratio of dust to HI 
gas. We adopt the scattering and absorption cross sections due to silicates and graphites derived 
by Draine & Lee (1984). The total extinction cross section is the sum of the contributions from 
scattering and absorption, and the relative proportion of these components is chosen so as to fit 
best the observed Galactic extinction law. Figure shows the sum a^xt of the extinction cross 



sections from graphites and silicates as a function of wavelength, assuming a dust-to-gas ratio of 
1 : 100 by mass (Whittet 1992). 




A. [/um] 



Figure 3: Sum of the extinction cross sections for scattering and absorption by graphites and 
silicates. The total cross section, croxt(A), is given in units of cm^ per 10^° hydrogen atoms, 
assuming a gas-to-dust fraction of 1 : 100 by mass, following Draine & Lee (1984). 



The dust optical depth is given by 



Text (A) = iVacxt(A) 



so that the extinction in magnitudes is 



Am(A) - |2.5 logio{exp[-7Vaext(A)]}| « 1.09 7Vaext(A) 



fll) 



(12) 



In quantifying the influence of dust, we restrict our attention to the observer's Johnson B 
band. In the rest frame of a lens at a redshift zi, this band is centered on the wavelength Ai = 
0.435 Aim (1 + zi)-\ For zi = 0.3, Figure | yields a^^tN w 0.084 iVao- Hence, for N = lO^i cm^^ 
Am « 0.9 magnitudes in the B band. Note that N can easily exceed this value for edge-on disks, 
which are favored both geometrically and due to the magnification bias. The value of N can be 
as large as iVo R/H = N^/q^ « 3.4 x 10^^ cm"^, yielding Am « 31 for a lens at zy = 0.3. This 
implies that the influence of dust on the imaging properties of spiral lenses can be severe. 

One would expect the dust-to-gas ratio of spiral disks to decline with increasing redshift, in 
accordance with their metallicity history. Indeed, Pei, Fall, & Bechtold (1991) infer that the dust- 
to-gas ratio in damped Lya absorbers at redshifts 2 < z < 3 is only a tenth of the Milky- Way 
value, although with considerable scatter (Fall & Pei 1993). If we parameterize the dust content 
to have a power-law dependence on redshift, (1 -\- z)^\ then a reduction by a factor of ^^ 10 from 
the present time to z w 2.3 implies 5 ~ 2. Equation (nj) is then changed to 



Am(A, zi) w 1.09 (1 + zi)-2 N20 CTcxt 



A 



z\ 



(13) 



This dependence will be referred to as the "dust evolution" model in the discussion that follows. 



2.4 Scalings of Galaxy Properties with Luminosity 

We assume that the lens population admits the Schechter luminosity function at present with a 
number density per unit luminosity, 



dnil) 



dl — n, l'^ exp(— ^) dl 



(14) 



where I = L/L^, is the scaled galaxy luminosity. For spiral galaxies, n^ — 1.5 x 10 "^h^Mpc '^ 
and v — —0.81 (Marzke et al. 1994). The luminosity is related to the circular velocity through the 
TuUy-Fisher (1977) relation, 

^ =?!/". (15) 

where a varies between ^ 2 in the B band and ~ 4 in the H band (see review by Strauss & Willick 
1995, and references therein). 

Equations (o) and (n5[) imply that the Einstein radius scales as & oc P/". The computational 
effort required later on is substantially reduced if the scale radii of the disk+halo model change 
with / in the same way as b does. For convenience, we therefore adopt the infrared value of a = 4 
and assume that the scale radii behave liken r^ oc Z^/^. It is then sufficient to calculate all cross 
sections for one reference luminosity only, and later scale the result to the desired I by changing 
the Einstein angle. If a^ is some arbitrary cross section for the reference luminosity i*, then the 
average cross section for the entire population of spiral galaxies is 
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2.5 Evolution of the Spiral Galaxy Population 

The simplest assumption about the evolution of the spiral population is that they maintain a 
constant comoving number density, n{z) = n(0)(l + z)'^, and constant scale radii. This is the 
no-evolution model. Alternatively, Mo, Mao, & White (1997) recently suggested a model for the 
evolution of spirals, based on four assumptions: (i) the disk mass is a fixed fraction m^ of the halo 
mass; (ii) the disk angular momentum is a fixed fraction jj of the halo's angular momentum; (iii) 
the radial disk profile is exponential, and the disk is ccntrifugally supported; and (iv) the disk is 
dynamically stable. These assumptions yield a set of simple and straightforward scaling relations, 
which we reproduce from Mo et al. (1997). 

We assume that the halo is a singular isothermal sphere with a radial density profile. 

Taking the virial radius r2oo as the size of the halo, the halo mass is M = v'^r2Qo/G. Here r2oo is 
the radius within which the average halo density is 200 times the critical density of the Universe, 

3M 3iJ2(z) , ^ 
— = 200 X ^-A^ , (18) 
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and hence, 



-"'^^^ -jowr '""""- (^20^) i^y • "'' 

where Hq = 100 /i km s^^ Mpc^^ is the Hubble constant. Similarly, 

M{z) = ^1— ^2.5x10^^ h-'M^f ^^ -Y f^W) . (20) 

^ ' lOGH(z) ° V220kms-i/ \ Hq J ^ ^ 



^Since r^j oc l^'^ implies constant surface brightness, the chosen scaUng reflects Freeman's law (Freeman 1970; 
Holmberg 1975; Peterson, Strom, & Strom 1979; van der Kruit 1987; Lauberts & Valentijn 1989). Thus, we 
efi'ectively ignore the scatter in Freeman's law, and the influence of low surface brightness (LSB) galaxies on the 
lensing cross section (e.g., Bothun, Impey, & McGaugh 1997 and references therein). This omission is justified here 
because (i) the average circular velocities of LSB galaxies is somewhat lower than that of high surface brightness 
(HSB) galaxies, and the lensing cross section is a sensitive function of the circular velocity, and (ii) LSB galaxies 
are less compact than HSB galaxies and thus less efficient lenses. For the purposes of strong lensing, it therefore 
seems safe to neglect any contribution from LSB galaxies. 



The disk radius r^ is related to the spin parameter of the halo, A (= J\E\^/^/GM^^'^, where J 
and E are the total angular momentum and energy of the halo, respectively), by 

and the central surface mass density of the disk isg 
10 mdVcH{z) f rrid^ 



So 






6.8 X 10" 



In summary, scale radii scale like oc VcH{z)~^, and the surface mass density^ scales like 
ex Vc H{z). Since the characteristic values of A, jd, and nid are expected to depend very weakly on 
redshift, we use these simple proportionality relations with Vc and H{z) in scaling the properties 
of the local spiral population to higher redshifts. Since H{z)/Hq > 1 for 2; > 0, scale radii decrease 
with increasing z, while the surface mass density increases with increasing z. 

The evolution in the number density of spiral galaxies can be expressed in terms of the Press- 
Schechter distribution function for the mass M, given by equation (pO[). We evolve the number 



density of galaxies given by equations (14), (|15|), and (20) with the factor nps(M, z)/nps(M, z — 
0), where nps(M, z) is the Press-Schechter number density of halos with mass M at redshift z. 
Because of the inherent uncertainty in modeling galaxy evolution, we present numerical results 
for three models, assuming: (i) no evolution; (ii) evolution of scale radii only; and (iii) evolution 
of galaxy number density and scale radii. 

It is well known that the present-day Press-Schechter mass function extends out to halo masses 
which are well beyond the galactic mass scale (e.g., Navarro, Frenk, & White 1995) in all viable 
models of structure formation. This is due to the fact that the nonlinear mass scale at present 
(which defines the exponential break in the Press-Schechter mass function) corresponds to a much 
larger circular velocity than L^ in the Schechter function does, based on the TuUy-Fisher relation 
[cf. Eqs. (|lj) and (|l^)]. This implies that non-gravitational processes (such as inefficient cooling, 
or expulsion of gas by supernovae) prevented disk formation inside super-galactic halos at the 
present time. An evolutionary model for galaxies based solely on the Press-Schechter approach is 
therefore incomplete. The simplest interpretation of the discrepancy between the Press-Schechter 
mass function and the local luminosity function of galaxies is that massive galaxies with L > L* 
did not change their dynamical properties since the redshift of galaxy formation (z ~ 2 — 4) when 
the nonlinear mass scale was comparable to their mass. Indeed, recent observations imply no 
significant evolution in the population of massive galaxies out to redshifts z ~ 1 (Ellis 1997, and 
references therein; but see Kauffmann, Chariot, & White 1996). However, the same observations 
reveal many more dwarf galaxies at high redshift than found locally. Since the lensing probability 
is dominated by L» galaxies at z < 0.7, the no-evolution model might be more appropriate 
for calculations of lensing by disk galaxies. However, to bracket the other extreme of complete 
evolution we show results also for the Press-Schechter prediction. Since we use the Press-Schechter 
mass function only to correct the overall normalization of the number density of spiral galaxies 
relative to its present-day value, our approach should be less affected by the incompleteness of the 
Press-Schechter treatment. 

In all numerical calculations, we use the cosmological parameters JIq = 0.3, f^A = 0, and 
h = 0.7. For the Press-Schechter mass function, we use the standard CDM power spectrum 

^Different realizations of galaxy formation in different tidal environments could result in disks of different sizes 
for the same halo properties, due to variations in the spin parameter acquired by the baryons. As discussed before, 
we ignore this scatter in the disk properties here. 

''We ignore the evolution of the HI mass fraction of the disk due to star formation, as most of the star formation 
occurred at redshifts 1.5 ^ ^i ^ 2.5 (Madau 1997), while most of the lensing probability is contributed between 
0.3<zi<0.7. ~ "^ 



(Bardeen et al. 1986) with the normahzation crgf^-i-j^p^ — I. A cosmological model with these 
parameters reproduces the local abundance of rich galaxy clusters (White, Efstathiou, & Frenk 
1993; Eke, Cole, & Frenk 1996; Viana & Liddle 1996), has the shape parameter Qoh ^ 0.21 
preferred by analyses of galaxy clustering (Peacock & Dodds 1994), but has a somewhat higher 
normalization than derived from the COBE data (e.g. Ratra et al. 1997). It also agrees with the 
observed abundance of giant luminous arcs in galaxy clusters (Bartelmann et al. 1997). 

3 Lensing Cross-Sections 

Next, we proceed to calculate the magnification cross section of the disk+halo lenses. This calcula- 
tion must be done numerically. First, we cover the source plane with a grid of source positions. Far 
from the lens center, the resolution of this grid can be low, while close to the caustic curves where 
the highest magnifications arise, the resolution should be high. We therefore use an adaptive grid 
in the source plane whose resolution increases towards the caustic curves. Then, for each source 
position, all image positions need to be found. For this purpose, we use the algorithm described 
by Schneider, Ehlers, & Falco (1992). Briefly, it is based on covering the image plane with an 
uniform grid, which is then mapped back to the source plane. All grid cells on the image plane 
whose mapping on the source plane contains the source position, are taken to contain images of 
the source. In regions of strong lensing, the parity of the cell might be flipped when it is mapped 
into the source plane, i.e. corners of the cell might be interchanged. One should therefore start 
with triangular rather than rectangular cells in the image plane, because the image of a triangle 
remains a convex figure whose interior is well-defined. This can be easily achieved by splitting 
each rectangular cell along one of its diagonals. Here again, we use an adaptive-grid approach in 
order to achieve high resolution at a reasonable computational cost. First, the image positions 
are searched on a coarse grid, constraining the regions on the lens plane where the images are 
located. Then, each of these regions is covered with a fine grid on which the final image positions 
are localized. This way, all the images corresponding to every segment on the source-plane grid 
are identified. 

The magnification of each image can then be computed from the total lensing potential (|9|), 



MGL = det-M^..-^-^), (23) 



evaluated at each of the image positions Xi. When dust extinction is included, the effective 
magnification of an image at position Xi is 

fi{xi) = ^iG-L{xi) e^p[-N{xi)acxt] , (24) 

where N{xi) is the HI column density at the position of the image. The net magnification is the 
sum of the moduli of the effective magnifications for all of its N images, 

N 

/. = ^|/i(f.)l- (25) 

This approach yields maps of the source magnification in the source plane. The magnification 
cross sections can be extracted from these maps as the area in the source plane within which 
sources are magnified by a factor > /x. The magnification cross section depends on the orientation 
of the disk relative to the line-of-sight. Denoting the cross section for a magnification > fj, and a 
disk inclination angle i by a'{fi,i), we obtain the inclination-averaged cross section through the 
integral, 

''■(m) — di cos{i)a'{fi,i) . (26) 

Jo 

Obviously, crip) depends on the cosmological distances involved and on the other lens parameters. 
We suppress these dependences here for brevity. 
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Figure 4: Inclination-averaged magnification cross section of a disk+halo lens (solid line) compared 
to that of a singular isothermal sphere (dashed line) with the same circular velocity. The source 
redshift is Zs = 2, and the lens redshift is zi = 0.3. The cross sections are given in units of the 
Einstein disk area, irb^. At small magnifications {n ^ 1) the cross section of the spiral lens becomes 
difficult to evaluate numerically, as indicated by the dotted section of the solid curve. The cross 
sections are almost identical, except for high magnifications, where the disk+halo model has a 
slightly larger cross section than the singular isothermal sphere. 

Figure m shows that the inclination-averaged magnification cross section of the disk+halo model 
is almost identical to that of a singular isothermal sphere with the same circular velocity. Without 
imposing further conditions, like a minimum image separation or a maximum flux ratio between 
the images, the lensing statistics of spiral lenses is well described by the simple spherical halo 
model. We will demonstrate below what happens when further constraints are imposed on the 
image properties. 

In the limit of small magnifications (/i ^ 1), the cross section for the disk+halo model becomes 
numerically incomplete; this follows from the fact that the cross-section increases rapidly as fi —^ 1, 
while the simulation is spatially bounded. This limit is irrelevant for our purposes because we 
always impose further imaging constraints. In particular, when we require the image separation 
and flux ratio to be bounded to reasonable limits, the cross section is confined to the multiple- 
imaging region, which is entirely contained within our simulated section of the source plane. 



3.1 Imaging Probabilities 

Given the inclination-averaged cross sections, the probability for a (point) source at redshift Zg 
to be imaged with magnification > /i is obtained by the line-of-sight integral over the density of 
lenses times their lensing cross section, cr(/i). For the model without number-density evolution, 
this integral yields 



^gl(m) 



■ bini 



A/a) 



dz(l + z)M§^ 



cdt 



dz 



ct(m) , 



(27) 



where |cdi/dz| is the proper-distance interval corresponding to the redshift interval dz. The factor 
n^blT{l + v + 4/a) comes from integrating over the luminosity distribution of the spiral galaxies 
[cf. Eq. (|l6|)]. When evolution of the galaxy population is taken into account as described in §2.5, 
the factor (1 + z)^ above is changed to nps(-M*, z)/nps(M*, 0) (1 + z)^, where M^,{z) is the mass 
of a galaxy with circular velocity Vc* at the corresponding redshift. 

Now let |dA^Qso/dS'|(5')d5' be the intrinsic number density of quasars at redshift Zs with a flux 
between S and S + dS. Accounting for magnification bias, the number of lensed quasars is 



%soiS) 



d^'^GL 



diV, 



QSO 



dS 



iS'), 



(28) 



and hence the probabihty for a quasar at redshift Zg to be detected with flux > S is 

1 



Pgi.{S) 



d5'PA.(4]^^(^') 



We approximate the observed quasar number counts in the i3-band by a broken power-law, 

diVqso _ r {S/S^Y for S < S^ 



AS 



(S/So)" for S> So ' 



(29) 



(30) 



where 5*0 corresponds to -Bqso ~ 19- For quasar redshifts Zg ^ 2, a — —1.64 and b = —3.52 (e.g. 
Pei 1995). 
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Figure 5: Probabihty ^^glCm) ^^r a quasar at a redshift Zg = 2 to be multiply imaged with 
magnification > /i, image separation > 0.3", and flux ratio of the images < 20. The curves 
are for the disk+halo model without dust (solid line), with non-evolving dust (dotted line), and 
with evolving dust (short-dashed line). We compare these results to a mass model of a singular 
isothermal sphere with no dust and the same asymptotic circular velocity as the disk-|-halo model 
(long-dashed line) . These curves do not include evolution of scale lengths or number density of the 
lens population. Here and in all the following figures, results for the singular isothermal sphere 
model are calculated without dust extinction. 

Figure |g shows PQL(/i), and Figure^ shows Pgl('5') for the no-evolution galaxy model. These 
figures include curves that illustrate the influence of non-evolving and evolving dust. All curves 
in Figure H were calculated under the additional constraints that the the two brightest images 
be separated by > 0.3", and that their flux ratio be < 20. The solid and long-dashed curves 
in Figure o are for the disk-|-halo model and for the singular isothermal sphere, respectively. 
Evidently, the disk-|-halo model is less efficient at producing multiple images with the specified 
properties than the singular isothermal sphere is, despite the fact that the total magnification 
cross sections of the two models are close to each other (Fig. 0). The reason for this difference is 
that the images produced by the spiral lens model are typically closer to each other than those 
of the singular isothermal sphere, and are occasionally below the 0.3" threshold. This occurs 
because a significant fraction of the lensing cross section is contributed by the disk in nearly edge- 
on orientations. The image separation on either side of the disk is then smaller than the Einstein 
diameter of the corresponding singular isothermal sphere. 

The imposed constraints on the minimum image separation and maximum flux ratio approxi- 
mately reflect characteristic thresholds for detecting a multiply imaged system through space-based 
observations with finite resolution and dynamical range (cf. Kochanek 1993). Figure ^ demon- 
strates that these additional constraints reduce the detection efficiency of spiral lenses considerably 
compared to the singular isothermal sphere model. The flgure also quantifles the severe effect that 
dust extinction has on the detection efficiency. Besides suppressing the detection probability by 
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Figure 6: Probability -Pgl(Sqso) for a quasar to be multiply imaged with image separation > 0.3" 
and flux ratio < 20, as a function of Sqso- The mass model includes no-evolution in panel (a), 
evolution of disk scale lengths in panel (b), and evolution of scale lengths and number density 
in panel (c). The four curves per panel are for the disk model without dust (solid curve), with 
non-evolving dust (dotted curve), and with evolving dust (short-dashed curve). The long-dashed 
curve shows results for the singular isothermal sphere model. 

one or two orders of magnitude for evolving and non-evolving dust, respectively, the dust extinction 
allows for values of fi — (S'obsorvcd/'S'intrinsic) Smaller than unity, which are otherwise unphysical. 

Figure g then shows the probability Pgl('S') for quasars at redshift Zg = 2 to be lensed with 
image separation > 0.3" and flux ratio < 20 as a function of the quasar B magnitude. Again, the 
dotted and short-dashed curves are calculated including dust extinction, and the solid and long- 
dashed curves are without dust for the disk+halo model and for the singular isothermal sphere 
model, respectively. The three panels in the flgurc are for different levels of cosmological evolution 
as detailed in the figure caption. The figure shows that bright quasars, with -Bqso < 19, are 
much more likely to be multiply imaged than faint quasars, due to the magnification bias. At 
-Bqso ~ 18, the multiple-imaging probability is about a factor of five higher than for -Bqso ^ 19. 
As already shown in Figure ||, dust extinction severely reduces the imaging probability, and it also 
leads to a weaker rise in the lens detection probability with increasing quasar brightness. 

3.2 Image Statistics 

In order to appreciate the significance of the selection effects in identifying spiral lenses, it is 
instructive to examine the distribution of image separations and fiux ratios. Figure ^ shows the 
cumulative distribution of multiple images as a function of the minimum image separation ^min- 
The distribution is arbitrarily normalized to unity at 6*111111 = 0.1". As before, the three panels in 
the figure are for disks with no-evolution (panel a), with evolution of scale radii only (panel b), 
and with evolution of scale radii and galaxy number density (panel c). The three curves per panel 
are for models without dust (solid line), with non-evolving dust (dotted), and with evolving dust 
(short-dashed). The long-dashed curve in each panel shows the singular isothermal sphere model 
without dust for comparison. 

Figure ^ implies that the median image separation is larger for the singular isothermal sphere 
than for the disk-|-halo model. Images produced by the spiral lens model are on average closer to 
each other by ~ 0.2" than expected from the singular isothermal sphere model. This is because 
multiple images in the dust-|-halo model are predominantly produced by the disk rather than by 
the halo, with a preference for nearly edge-on disk orientations. 

Dust increases the average image separation because images with small separations occur close 
to the disk, where their flux is heavily suppressed by extinction. Disk evolution reduces the average 
image separation slightly, because the disks get more compact with increasing redshift. Number 
density evolution according to the Press-Schechter model increases the abundance of fleld spirals 
at redshifts z ^ 1 — 2 (since some of these galaxies have subsequently merged and disappeared by 
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Figure 7: Cumulative distribution of multiple images with separation 6 > 0,nin, as a function of 
^min, for a constant maximum flux ratio rmax = 20 and quasar magnitude iJqso = 18- The disk 
model includes no-evolution in panel (a), evolution of scale lengths in panel (b), and evolution of 
scale lengths and disk number density in panel (c). The three curves per panel are without dust 
extinction (solid line), with extinction by non-evolving (dotted) and evolving (short-dashed) dust. 
The long-dashed curve describes the singular isothermal sphere model with no dust. 

now), and reduces the image separation further because the lenses are on average further away. 
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Figure 8: Cumulative distribution of multiple images with a B band flux ratio r < 



<;, as a 



function of r 

same meaning as in the previous figure. 



for a constant 6*1111,1 = 0.3" and Bqso = 18. The panels and the curves have the 



Figure S shows the cumulative distribution of images as a function of the maximum flux ratio 
?'max, arbitrarily normalized to unity at rmax = 50. The notation is the same as in Figure R. Gen- 
erally, the curves for the disk-f-halo model without dust are the flattest, showing that the images 
are usually of comparable brightness. When dust is included, the average flux ratio increases due 
to the strong extinction gradient around the disk. 

4 Effects of the Gaseous Disk 

4.1 Impact of Dust on Lensing Statistics 

To further examine the influence of dust on the detectability of lensed quasars, we show in Figure S 
the ratio between the lensing probability -Pgl(Sqso) with and without dust. The three panels 
are again for a galaxy model with no evolution (panel a), with length scale evolution only (panel 
lb), and with length-scale plus number-density evolution (panel c). The two curves per panel are 
for non-evolving and evolving dust (dotted and solid lines, respectively). 
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Figure 9: Reduction in the lens detection probability due to dust extinction in the (observer's) B 
band. We assume 6*111111 = 0.3" and r^^x = 20. The curves show the lensing probabiHty including 
dust divided by the lensing probability without dust, as functions of the quasars' B magnitude. 
The disk model has no-evolution in panel (a), evolution of scale lengths in panel (b), and evolution 
of scale lengths and disk number density in panel (c) . The two curves per panel are for non-evolving 
(dotted line) and evolving (solid line) dust. 

Figure ^ shows that the majority of all quasars which are lensed by spirals are undetectable 
in the B band. The deficit caused by dust increases for brighter quasars as a result of the 
magnification bias. At Bqso ~ 18, only about 10% of all lensed quasars are observable, quite 
independent of the degree of evolution that is included in the calculation. If dust does not evolve 
with redshift, the deficit of lensed quasars is more severe. For faint quasars, Bqso ^ 20, the 
reduction in the lens detection probability is least severe in the no-evolution model. This is 
mainly because when disk evolution is included, the surface density of the disk increases with z, 
leading to more extensive dust extinction than the no-evolution model predicts. 

Figure ^ implies that the fraction of quasars lensed by spiral galaxies in radio surveys should 
be higher by about an order of magnitude relative to that found in optical surveys. 

4.2 Statistics of Damped Lya Absorption 

Based on the neutral hydrogen column density at the position of each image, we can investigate 
the influence of lensing and dust on the statistics of damped Lyman-a absorption by spirals. Let 
Pgl(»S', N) be the probability to observe a quasar with a flux > S which is imaged by the population 
of spirals and shows an HI column density > TV in its spectrum. When there is more than one 
image, we take N to be the HI column density in the brightest image, i.e. that which dominates 
the absorption trough. As shown by Bartelmann & Loeb (1996), the observed column-density 
distribution of neutral hydrogen is then given by 



^(^) = i^ 



dPGL{S,N) 



ON 



(31) 



where AX is the absorption distance scanned by the damped Lya absorbers (DLAs) in the survey 
(Bahcall & Peebles 1969). Examples for the distribution of Nf{N) are plotted in Figue 11^. 
The inferred cosmological density parameter in neutral hydrogen, fiHij is given by 



r^HT — — 



Hq m 



C Pcr.O Jo 



dNN f{N) 



(32) 



where Pcr,o is the present-day critical density of the Universe, and m is the mean molecular mass. 
Because of the magnification bias, f{N) depends on the quasar magnitude Bqso, and so does the 
inferred rini- Figure O shows the ratio between the inferred and true values of r^Hi for spiral 
galaxies. The different panels and line- types are the same as those in previous figures. Similarly 
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Figure 10: The HI column-density distribution, Nf{N), for lensed quasars with _B = 18 at ^g = 2 
as a function of A^20 = (-/V/10^° cm"^). The disk model includes no-evolution in panel (a), 
evolution of scale lengths in panel (b) , and evolution of scale lengths and number density in panel 
(c). The three curves were calculated without dust, with non-evolving dust, and with evolving 
dust (sohd, dotted, and short-dashed curves, respectively). For reference, the long-dashed curve 
shows the result without lensing or dust extinction. 



to Pql, f{N) and f^Hi depend on the redshift range of the absorbers [cf. Eq. (|27|)]. We use the 
full redshift range, < z < Zg = 2 for Figure |l^ and for the upper panels of Figure |ll|, and the 
high-redshift range 1 < z < Zg for the lower panels of Figure ^ 

Models without dust cause Ohi to be overestimated by factors of up to a few in samples of 
bright quasars. The magnification bias brings into view quasars that were otherwise too faint to 
be detected. The lensed quasar images occur close to the lens center, where the HI column density 
is higher than average. However, when dust is included, the net effect is reversed. Now those 
quasars whose images are close to the disk are obscured, and the number of quasar spectra with 
high HI column density is lower. For evolving dust and Sqso ~ 18, only about 20%-30% of the 
neutral hydrogen is detected. 



5 Conclusions 

We investigated the lensing effects of spiral galaxies, modelled as maximal truncated Mestel disks 
embedded in spherical halos. This model makes the disk as massive as possible. Our results 
therefore represent the opposite extreme to the simple isothermal sphere models. Intermediate 
models with less massive disks, would shift these results towards those obtained from the spherical 
model. 

Although a disk-|-halo configuration yields an inclination-averaged lensing cross section similar 
to that of the spherically symmetric mass distribution (Fig. ^, its typical image separations are 
considerably smaller (Fig. 0). A substantial contribution to this cross section comes from lensing by 
nearly edge-on configurations of the disk, in which a pair of images straddles the disk. Because the 
images are often close to the disk, they suffer strong extinction by dust. The extinction lowers the 
detection efficiency of spiral lenses in the optical band by an order of magnitude. Moreover, disk 
lenses are often characterized by small image separation (Fig. 0) and large differential extinction 
(Fig. g). As a result, the selection effects imposed by the finite angular resolution and dynamic 
range of observations set limits on the minimum image separation and their maximum flux ratio, 
and lower substantially the probability for observing spiral lenses (cf. Fig. ||). 

The spiral-lens system B 0218+357 shows compelling evidence for strong extinction. O'Dea 
et al. (1992) find that the quasar spectrum is red; Wiklind & Combes (1995) and Menten & 
Reid (1996) find strong molecular lines; and Grundahl & Hjorth (1995) find that image A is much 
fainter than image B in the optical in contrast to the radio observations, arguing for a substantial 
extinction of image A. Jaunsen & Hjorth (1997) have recently argued for the existence of disk 
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Figure 11: The ratio between the inferred and the true values of the density parameter of neutral 
hydrogen as a function of the quasar B magnitude. For the upper three panels, the redshift range 
for the absorbers is < z < Zg = 2, while for the lower panels it is restricted to 1 < z < Zg. The 
disk model includes no-evolution in panels (a) and (d), evolution of scale lengths in panels (b) 
and (e), and evolution of scale lengths and number density in panels (c) and (f). The three curves 
were calculated without dust, with non-evolving dust, and with evolving dust (solid, dotted, and 
dashed curves, respectively). 

lensing and associated extinction also in the lens system B 1600-1-434. More generally, Malhotra, 
Rhoads, & Turner (1997) argue for a systematic reddening of lensed quasars relative to the rest 
of the quasar population. 

The optical depth to gravitational lensing of quasars was shown to be a very effective tool 
for setting constraints on the cosmological constant (Kochanek 1996, and references therein). In 
singular isothermal models, ^ 80% — 90% of the lensing probability is contributed by elliptical 
galaxies. However, if a significant fraction of the present-day ellipticals were in the form of spiral 
building blocks at z ^ 1 (as argued by Kauffmann, Chariot, & White 1996), then the depletion 
of spiral lenses due to dust extinction and selection effects could weaken the current lensing 
constraints on the cosmological constant. This issue was addressed by Mao & Kochanek (1994) 
and Rix et al. (1994). They found that estimates of the cosmological constant from the statistics 
of strong lensing can be significantly changed only when the elliptical galaxies seen today evolved 
dramatically below a redshift z ^ 1. 

In this context, and in view of ongoing lensing surveys, it is instructive to examine which 
fraction of the total lensing optical depth is contributed by spirals rather than by ellipticals. To 
calculate this, we assume that the total galaxy population today is composed of 25% ellipticals and 
75% spirals. We model spirals as described above, and ellipticals as singular isothermal spheres 
with a characteristic velocity dispersion of CTi,* = Wc*/-\/2 = 220kms~^, following the Schechter 
luminosity function (HJ) and the Faber-Jackson (1976) relation. For simplicity, we use the same 
Schechter-function parameters as for the spirals. We consider two scenarios, one in which the 
number density of ellipticals stays constant with redshift, and another in which ellipticals are 
assembled by merging of spirals. In the latter scenario, we assume that the number density of 
ellipticals changes with redshift as a power-law, n^(z) — nE(0) (1 -I- zf ^ with e w —1.6, chosen 
such that ^ 2/3 of the ellipticals are in the form of spiral bulding blocks at z ~ 1 (Kauffmann 
et al. 1996). We further assume that mergers conserve mass. In both scenarios, we can then 
compute the lensing optical depths for both spirals and ellipticals ^^ql-e s fron^ equation ( |27| ) and 
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the corresponding imaging probability Pgl:E,s from equation (E9f) (under the constraints that the 
images be separated by > 0.3" and have a flux ratio < 20). 

Figure R2 shows the fraction of the total lensing probability contributed by spirals, namely 
(■PGL,s/^GL,totai) = ^gl,s/(-Pgl,e + ^GL,s)- Panel (a) shows resuhs for eUipticals with constant 
comoving number density, while panel (b) examines the scenario in which ellipticals merge out 
of spirals. The three curves per panel are for spirals without dust, with non-evolving dust, and 
with evolving dust (solid, dotted, and dashed lines, respectively). The no-dust results are valid for 
radio-selected lenses. The solid line in panel (a) reproduces the familiar result that dust-free spirals 
contribute at most ~ 20% of the multiple-imaging probability, even for the small image-separation 
cutoff of 6'min = 0.3"; this fraction is lower for higher 6'niin- In the presence of dust, the fraction 
contributed by spirals is considerably smaller. In the merger scenario, the spiral lens fraction rises 
to ~ 50% without dust, and '^ 9% with evolving dust, for quasars with B ~ 18. Thus, up to a 
half of the lenses found in radio surveys (which are sensitive to small image separations) might be 
imaged by spirals. 
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Figure 12: Fraction of the total multiple-imaging probability contributed by spiral galaxies. In 
panel (a) the elliptical population is assumed to have constant comoving number density, and 
in panel (b) the number density of ellipticals evolves as a power law of redshift so that 2/3 of 
the ellipticals are in the form of spiral bulding blocks at z = 1. The three curves per panel are 
for spirals without dust (solid line), with non-evolving dust (dotted line), and with evolving dust 
(dashed line). Apart from merging to ellipticals, any other evolution of the spiral population was 
ignored. 

Figure ^ implies that radio surveys should be about a factor of 5-10 more efficient at detecting 
spiral lenses than optical surveys. The recent CLASS and JVAS radio surveys indeed provide 
preliminary hints of a more substantial population of spiral lenses than found in optical surveys 
(cf. Table 1 in Browne et al. 1997; and Jackson, Nair, & Browne 1997; but see Fassnacht & Cohen 
1997). To date, the combined CLASS and JVAS samples encompass more than lO'^ flat-spectrum 
radio sources, 11 of which have been identified as lens systems. A substantial fraction of those 
has been classified as being lensed by spirals or SO galaxies (cf. Myers 1997; Williams & Schechter 
1997). Because of their small image separation and simple geometry, lenses containing edge-on 
disks might offer a unique opportunity for estimating the masses of galactic disks at high redshifts, 
and also for constraining the Hubble constant based on the time-delay between the flux variations 
in their images. Indeed, the time delay in B 0218 has been measured to 12 ±3 days, translating to 
a Hubble constant of Hq ~ 60kms~^ Mpc~^ (Corbett et al. 1996), and B 1608+656 offers another 
promising lens system for this purpose (Myers et al. 1995; Fassnacht et al. 1996). 

Finally, we find that dust dominates over the magnification bias due to lensing in modifying 
the statistics of damped Lya absorption by spirals (Fig. O) . If most of the damped absorption 
systems are spirals, the inferred value of the cosmological density parameter in HI at z < 1 could 
be underestimated by a factor of a few. 

The above results were obtained under the assumption that all spirals possess the same dust- 
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to-gas ratio throughout their disks at any given redshift. If the scatter in the dust content of 
different galaxies is large or if the dust distribution is patchy, then the quantitative impact of dust 
on the lensing statistics would be altered. In addition, if spiral disks at 2; < 1 are less massive 
than our maximal disk model assumes, then the lensing cross section of the disk and its related 
extinction signature would be reduced. 
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